Preprint typeset in JHEP style - HYPER VERSION ULB-TH/09-21 



IDM & iDM 
or 

The Inert Doublet Model and Inelastic Dark Matter 



Chiara Arina, Fu-Sin Ling, Michel H.G. Tytgat 

Service de Physique Theorique, Universite Libre de Bruxelles, CP225, Bid du Triomphe, 
1050 Brussels, Belgium 

E-mail: carina@ulb.ac.be, fling@ulb.ac.be, mtytgat@ulb.ac.be 

Abstract: The annual modulation observed by DAMA/Nal and DAMA/Libra may be 
interpreted in terms of elastic or inelastic scattering of dark matter particles. In this paper 
we confront these two scenarios within the framework of a very simple extension of the 
Standard Model, the Inert Doublet Model (IDM). In this model the dark matter candidate 
is a scalar, the lightest component of an extra Higgs doublet. We first revisit the case for 
the elastic scattering of a light scalar WIMP, Mdm ~ 10 GeV, a scenario which requires 
that a fraction of events in DAMA are channelled. Second we consider the possibility of 
inelastic Dark Matter (iDM). This option is technically natural in the IDM, in the sense 
that the mass splitting between the lightest and next-to-lightest neutral scalars may be 
protected by a Peccei-Quinn (PQ) symmetry. We show that candidates with a mass Mdm 
between ~ 535 GeV and ~ 50 TeV may reproduce the DAMA data and have a cosmic 
abundance in agreement with WMAP. This range may be extended to candidates as light 
as ~ 50 GeV if we exploit the possibility that the approximate PQ symmetry is effectively 
conserved and that a primordial asymmetry in the dark sector may survive until freeze-out. 
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1. Introduction 

DAMA/Libra, and the former DAMA/Nal are dark matter (DM) direct detection experi- 
ments (DAMA in the sequel) which have observed 11 successive cycles of annual modulation 
in the rate of nuclear recoils, with a statistical significance of 8.2a (1). These measure- 
ments are consistent with the signal that would arise from collisions between nuclei and 
weakly interacting massive particles (WIMP) from a galactic halo, the flux of dark matter 
particles being modulated by the periodic motion of the Earth around the Sun (2; 3). For 
a standard, Maxwellian, spherical, non-rotating halo, the collision rate is expected to peak 
around the 2nd of June, when the velocity of the Earth with respect to the halo of dark 
matter adds to the velocity of the Sun in the Galaxy. 

A problem with DAMA is that a straightforward interpretation of the data in terms 
of WIMP collisions is in contradiction with the exclusion limits set by other direct de- 
tection experiments (see for instance (4) or (5) for recent reviews). Ways to reconcile 
the DAMA observations with the other experiments include models of mirror dark matter 
(6; 7), composite dark matter (8), MeV dark matter (9), hidden charged dark matter (10), 
dipolar dark matter (11), etc. In the present work, we focus on two alternative scenar- 
ios. Concretely, we consider the spin-independent (SI) elastic scattering of a light WIMP, 
with My)m in the multi-GeV range (12-25), and the so-called Inelastic Dark Matter (iDM) 
scenario (26-31). To a large extent, our conclusions concur with those obtained in these 
works, but our main intention is distinct, as we wish to confront a specific particle physics 
model, the so-called Inert Doublet Model (or IDM), to the recent DAMA data. 
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The IDM is a very simple extension of the Standard Model, first introduced in (32) 
but which has received some attention more recently, starting in (33; 34). Despite its 
simplicity, the IDM has an interesting phenomenology, which may be related to the origin 
of neutrino masses (34) or to electroweak symmetry breaking (33; 35) and various aspects 
of dark matter have been discussed (24; 36-41). Also, the IDM is one instance of a whole 
category of models for dark matter, to which models like Minimal Dark Matter (42) or 
Scalar Multiplet Models (43) belong. 

The IDM contains two Higgs doublets and, to prevent FCNC, a discrete Z<i symmetry 
is imposed, with the extra (or inert) Higgs taken to be odd under Zi. All the other fields, 
with the possible exception of right-handed neutrinos (34), are even. There are two extra 
neutral scalars states, noted Hq and Aq, which, if Z2 is unbroken, may play the role of dark 
matter. Here the lightest state, hence the dark matter candidate, is taken to be Hq. The 
mass splitting between these two states is controlled by a quartic coupling to the Higgs, 
here A5, with 

M 2 Ho - M 2 Ao = X 5 v 2 

where v = 246 GeV is the vacuum expectation value of the Higgs. In the limit A5 — ► 0, the 
Z2 symmetry is elevated to a global U{\)pq Peccei-Quinn symmetry, a limit relevant for 
the inelastic dark matter scenario, which requires a small splitting, Mh — Mj\ ~ 100 keV. 
In this limit, scattering and annihilation are dominantly through the Z channel. Another 
interesting limit is that of a light Hq and two heavier, nearly degenerate states, i.e. Aq and 
the charged partner in the inert doublet. This limit, which is relevant for elastic scattering 
of a light WIMP scenario, is also protected by a global symmetry, which is SU(2) in this 
case, analog of the custodial symmetry of the Standard Model. In this scenario, scattering 
and annihilation are dominantly through the Higgs channel. The elastic and inelastic 
scenarios for DAMA correspond to specific limits in the parameter space of the IDM. The 
main purpose of this paper is to discuss the behavior of the model in these two specific 
regions, in connection with the present direct detection experiments. 

The organization of the paper is as follows. The IDM is introduced in Section 2. The 
physics of direct detection, including the seasonal modulation, and the relevant data are 
discussed in Section 3. We discuss the elastic scenario in Section 4. The main result of 
this section is summarized in Figure 3, in which we compare the goodness of fit of a light 
WIMP candidate in IDM to the DAMA data, together with the exclusion limits set by 
the other experiments and with the WMAP relic abundance. At 3a, only a small region 
of the parameter space is compatible with all observations. We then discuss the inelastic 
scenario in Section 5. The main result of this section is summarized in Figure 8, in which 
we show that a whole range of candidates between Mdm ~ 535 GeV and Mdm ~ 50 
TeV is compatible with all observations, including the WMAP abundance. In the same 
section, we argue that lower mass candidates, which in the standard freeze-out scenario 
have a too small relic abundance, may be consistent with all observations, provided there 
is a primordial PQ asymmetry in the dark sector. We give our conclusions in Section 6. 
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2. The Inert Doublet Model 

The Inert Doublet Model (IDM) is a simple extension of the Standard Model with dark 
matter, with two Higgs doublets and a Z2 symmetry. The usual Brout-Englert-Higgs 
doublet (hereafter the Higgs doublet) is denoted by H\. The extra, or inert doublet, H2, is 
the only field of the model that is odd under the Z 2 symmetry. 1 This ensures the stability of 
the lightest member of H2, which will be a DM candidate, and prevents from flavor changing 
neutral currents (FCNC) (32). We will assume that Z2 is not spontaneously broken and that 
H2 does not develop a vacuum expectation value. We write H2 = (H + (Hq + iAq) / \/2) T , 
similarly to the ordinary Higgs doublet, and Hi = (h + (vo+h+iGo) / \/2) T . The potential 
is written as 

V(H U H 2 ) = rf\Hi\ 2 + v 2 2 \H 2 \ 2 + Ai|iii| 4 + A 2 |# 2 | 4 
+ X 3 \Hi\ 2 \H 2 \ 2 + X a \h\h 2 \ 2 + ^ 

After electroweak symmetry breaking, (Hi) = v = — /uf/Ai = 246 GeV, the masses of the 
physical scalar fields are given by 

M 2 h = 2X x v 2 , 

M H = f4 + ^H V 2 , 

Ml = nl + \a v 2 , 

M 2 H+ = li\ + X Hc v 2 , (2.2) 

with \h c = A3/2 and Xh ,Ao = (A3 + A4 ± As)/2. We will consider Hq to be the DM 
candidate (i.e. A5 < 0) though the results would be exactly the same for Aq changing the 
sign of A5. Notice that most of the negative couplings parameter space is excluded by 
vacuum stability constraints. Indeed, to ensure that the scalar potential is bounded from 
below, we need (33) 

Ai, 2 > , 

Aho , Aa , Xh c > -\AiA 2 . (2.3) 

The IDM has already been extensively studied in the literature. It has been shown 
that a viable DM candidate with the correct relic abundance can be obtained in three 
regimes, low-mass (M# <C mw) (24; 35), middle-mass (Mh mvi/)(33; 36) and high- 
mass (Mh 3> mv^)(36; 43). Direct and indirect detection constraints were investigated in 
Refs. (24; 33; 36-40; 46) and confrontation to colliders data and related future prospects 
was done in (47; 48). Here we will focus on the implications of DAMA for the IDM. 

Two particular limits of the mass relations will be most relevant here. The most 
obvious one is the limit in which A5 — > 0, in which case the neutral particles are degenerate, 

1 Actually, this is not the only logical possibility. Since fermions come in pairs, as one may also impose 
that all the fermions of the Standard Model are odd under Z2, without jeopardising the stability of the 
neutral scalar (44; 45). As emphasized in these works, this so-called 'matter-parity' naturally points to a 
scalar dark matter candidate and furthermore opens the possibility to embed the model in SO (10). 



{h\h 2 ) 2 + h.c. 
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Mh = Ma , and the theory is invariant under a larger, Peccei-Quinn symmetry, U(1)pq D 
Z2- This limit will be relevant in Section 5 where we will discuss the inelastic scenario 
for DAMA. The existence of an enhanced symmetry will imply that the limit of nearly 
degenerate neutral scalars, albeit fine-tuned, is technically natural. From Eqs. (2.2), a 
small mass splitting 5 = Ma — Mjj a — 100 keV as suggested by DAMA corresponds to a 
coupling 

^ ^(m^Xm) ■ <"> 

It is worth noticing that the IDM is the only model among Scalar Multiplet Dark Matter 
models, Ref. (43), where such a tiny mass splitting between the neutral components can be 
generated at the renormalizable level. This opens the interesting possibility that for higher 
scalar multiplets the splitting may be generated by higher order operators. 

Another, perhaps less obvious limit is that with Ma = M H ±, in which case three 
real states are degenerate and the theory is invariant under a global SU{2) symmetry 
(35; 49). This symmetry is analogous to the custodial symmetry in the Standard Model, 
in the sense that it prevents the appearance of large radiative corrections to the p, or 
equivalently T, parameter. This implies that, to some extent, it is possible to decouple 
the extra components of the inert doublet, Mh <C Ma ~ Mjj±. In this limit, the dark 
matter phenomenology is similar to that of a singlet scalar field (24), a model for dark 
matter much discussed in the literature (50-52). 



3. Direct detection experiments 

For the sake of reference, we give all the relevant steps necessary to compute the signal 
due to SI scattering of dark matter on nuclei in direct detection experiments. A most 
useful reference is the review by Lewin and Smith (53). See also Jungman, Griest and 
Kamionkowski (54). 



3.1 Event rate 

We consider collisions between dark matter particles from the Galactic halo and the nuclei 
of a given low background detector. The relevant characteristics of the halo are the local 
density of dark matter, taken to have the fiducial value pdm = 0.3 GeV/cm 3 at the Sun's 
location, and the local distribution of velocities f(t7) with respect to the Earth. 

At a given recoil energy En, the differential event rate of nuclear recoils can be factored 

as 

dR Pdm da 

dE~ R = M^dEn^t) (3.1) 

where Mdm is the WIMP mass, da/dE^ encodes all the particle and nuclear physics 
factors, and i](Er, t) is the mean inverse velocity of the incoming particles that can deposit 
a given recoil energy Er. The time dependence of the velocity distribution is induced by 
the motion of the Earth around the Sun, which leads to a seasonal modulation of the event 
rate (2; 3). 
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The total recoil rate per unit detector mass in a given energy bin [E\, E2] is obtained 
by integrating Eq. (3.1), 



F\E R ) , (3.4) 



R(t) = J* dE e(E) * G(E, a(E))^j , (3.2) 

where e is the efficiency of the detector and the finite energy resolution of the experiment 
is taken into account by convoluting the differential rate with a gaussian distribution with 
spread cr{E). For detectors made of several elements, the total rate is the average of the 
rates Ri{t) for each component i, weighted by its mass fraction fi 

R(t) = J2fiMt) (3-3) 

i 

Finally the expected number of observed events per unit time is the product of the total 
rate times the detector mass M^et- We will discuss further effects, like quenching factors 
and channelling effects in the next subsections, together with the other characteristics of a 
specific experiment. 

The particle physics is enclosed in the term da/dEn, which is generally parametrized 

as 

da = MxaZ (fv Z + ( A ~ Z )f n] 
dE R ~ 2^ /2 

where is the nucleus mass, \i n is the reduced WIMP/neutron mass, <r° is the zero 
momentum WIMP-neutron effective cross-section, Z and A are respectively the number of 
protons and the atomic number of the element, and f p (f n ) are the WIMP effective coherent 
coupling to the proton (resp. neutron). The nuclei form factor F 2 (Er) characterizes the 
loss of coherence for non zero momentum transfer. We use the simple parametrization 
given by Helm (53; 55), defined as 

F(E R ) = 3e~^ ^ 'J^ ^ (3.5) 

with q = \/2MnEr, s = 0.9 fm and r the effective nuclear radius. This form factor is 
optimal for scattering on Iodine (53). For simplicity we use the same form factor for all 
targets (more accurate form factor are off by at most 0(20%), for some targets and at 
large recoil energies (56)). 

Finally, the velocity distribution appears in the quantity 

V(E R ,t)= [ d 3 v^ , (3.6) 



where v the WIMP velocity wrt the Earth and v m %n is the minimum velocity needed to 
provoke a recoil inside the detector. The threshold velocity is given by 
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with // the nucleus- WIMP reduced mass. This formula encompasses both elastic (5 = 0) 
and inelastic (S 7^ 0) scatterings. Here 5 is the mass splitting between the two DM particles 
involved in the scattering (which are different in the case of inelastic scatterings). 

When the velocity distribution in the galactic frame f ga i(^) is isotropic, rj is given by 

rj = f + (F( Vesc ) - F(v)) dv , (3.8) 

with v± = min{t> esc , v m i n rb i>©}, v esc is the galactic escape velocity, u©(i) = vq + VEo(t) 
is the Earth's velocity in the galactic frame, with v Q = (0, 220, 0) + (10, 13, 7) km/s is the 
velocity of the Sun wrt the halo (57), t>Eo = 29.8 km/s is the Earth mean orbital velocity, 
and F(v) = J vf ggb \(v) dv. Notice that F(v) is an even function of v because i ga i(v) only 
depends on the modulus of the velocity. In deriving this formula, we have used the realistic 
assumption v esc > at any time t. For a standard Maxwellian distribution, 

fgal(^ = T77 Vlm^ e ~ v2,Vl ^iv<v esc , (3.9) 

where N(v esc ) is a normalization factor, Eq. (3.8) reads 

V = ATf \ (Erf ( V -±) - Erf (?=.) - 2( V+ ~ V ~ ] e-^'A , (3.10) 

2N(v esc )V (B { \V / \V J V7TU0 J 

which agrees with the result of Ref. (21). It might seem curious that Eq. (3.8) is expressed 
as a definite integral between v- and v+ while particles with any velocity above v m in should 
contribute to r/. However, the function F(v) is itself an integral on the velocity distribution. 
Also the dependence in v esc is explicit. 

In this paper, all the fits and x 2 analyses have been performed with a standard 
Maxwellian halo velocity distribution with a dispersion parameter vq = 220 km/s, and 
an escape velocity v esc = 450 and v esc = 600 km/s or v esc = 650 km/s, somewhat on 
the edges of the expected range, 498 km/s < v esc < 608 km/s (58). Modifications of the 
velocity distributions and their impact on the fits, have been discussed more in depth in 
various works, see for instance Refs. (20; 22; 29; 59). 

3.2 DAMA 

The former DAMA/Nal (60) and present DAMA/LIBRA (1) experiments are made of 
Nal(Tl) crystals. They are designed to detect the dark matter recoil off nuclei through the 
model independent annual modulation signature due to the motion of the Earth around 
the Sun. The experimental results obtained by DAMA/LIBRA, with an exposure of 0.53 
tonxyr collected over 4 annual cycles, combined with the ones of DAMA/Nal, for an 
exposure of 0.29 tonxyr collected over 7 annual cycles, corresponding to a total exposure 
of 0.82 tonxyr, show a modulated signal with a confidence level of 8.2 a (1). 

For an isotropic velocity distribution, the contribution of the signal to the counting 
rate in each bin can be approximated as 

dR 

— ~S (E R ) + S m (E R )coBu(t-to) (3.11) 
dbj R 
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where So is the average signal, S m is the modulation amplitude, u = 2-ir/T with T = 1 
year, and to is the phase corresponding to to = 152.5 day (June 2), the time of year at 
which t>e(t) is at its maximum. The amplitude of the modulated signal is therefore given 
by: 

„ , , dR 1 r dR , , dR ,^ , "i . 



DAM A gives the modulated average amplitude over bin intervals, with 

S m = p 1 jp dE S m (E) . (3.13) 
— Hj\ J Ex 

The energy resolution of the detector given by the collaboration (1) is 

a(E) 0.45 



E JE 



+ 0.0091 (3.14) 



and the efficiency is e = 1 (1; 21). 

The DAMA spectrum is given in keVee (electron equivalent eV). The observed en- 
ergy released in scintillation light is related to the nuclei recoil energy through a so-called 
quenching factor q, E sc ^ = q • E Teco ^. This expresses the fact that a recoiling nucleus 
may loose energy by collisions with other nuclei, hence in the form of heat, or through 
collisions with electrons, which create scintillation light. The reference values for Iodine 
and Sodium are respectively qi = 0.09 and qjy a = 0.3. However it has been pointed that so- 
called channelled events may play a role (17; 61). This refers to events in which a recoiled 
nucleus moves along the axis of the Nal crystal, losing most of its energy by collisions with 
electrons, in which case the quenching factor may be larger, up to q ~ 1. Once channelling 
is taken into account, collisions of light WIMPs with Iodine become relevant, while recoils 
on sodium are negligible in all scenarios (17-21; 62). We use the fraction / of channelled 
events given in Ref. (20), 

e -E R /18 e -E R /40 

ha{ER) = 1 + 0.75E R ' h = 1 + 0.6525* (3 - 15) 

where the recoil energy En is in keV. We have verified that the other parametrizations 
that we have found in the literature, (18; 21; 28), give identical results. 

We compare the prediction of the IDM to the observed energy spectrum of Eq. (3.12), 
using all the data provided by DAMA, which are given in Figure 9 of Ref. (1), in 36 bins 
of width 0.5 keVee, from 2.0 to 20.0 keVee (1). For the explicit numerical values, we refer 
to the Table III of Ref. (21). 

We use a basic statistical analysis. First we keep all the 36 bins, even though most of 
the signal is supposed to be concentrated around the first twelve to eighteen bins. This 
leads to some dilution of the statistical significance of the signal, which in turns means that 
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we are more optimistic about which parameters of our model may actually fit the data. 
However our results are essentially consistent with those of (19-21; 28; 29) 2 . 



We use a Goodness-of-Fit (GOF) method with a \ 2 metric to determine the allowed 
regions at 90%, 99% and 99.9% of confidence level (CL), with 



where Si are the theoretical predictions and S° the reported signals in each bin and <Ji 
are the experimental uncertainties in the measurements. For 36 bins and a scan over two 
parameters the 90% (99% and 99.9%) CL corresponds to x 2 < 45 (resp. x 2 < 56 and 



We have made numerous comparisons with the other existing analyses, and have found 
good agreement (19-21; 28; 29). As a further check, we have also repeated the above 
procedure considering only a two bins version of the data (i. e. the actual numbers published 
by DAMA) for 2-6 keVee and 6-14 keVee bins, as in Refs. (18; 20; 21). As is well-known, 
the allowed parameter regions are much larger but in the present work, so as to avoid the 
cluttering of the figures, we only show our results using the complete spectral data. 

As emphasized in Refs. (18; 20; 21), the total unmodulated rate provides a further 
constraint on the DM parameters, we have compared the total experimental rate from (1) 
with the unmodulated spectrum in each bin, and required that the predicted total number 
of events do not exceed the measured signal. 

3.3 Exclusion limits 

So far all the other direct detection experiments searching for dark matter are compatible 
with null results. In this section we briefly describe the experiments that lead to the most 
constraining limits on both the elastic, Section 4, and the inelastic scenarios, Section 5. 

Light nuclei, like Aluminium and Silicon, are more sensitive to light WIMPs scattering 
Mdm ~ multi-GeV, and provide the strongest upper bounds on the allowed parameter 
space favoured by DAMA. In the inelastic scenario, which involves heavier candidates, ex- 
periments made of heavy nuclei, like Iodine, Xenon and Tungsten are the most constraining 
ones. Germanium made detectors fall in-between. 

In computing the rate of Eq. (3.2) for each experiments, we uniformly assume that the 
small number of events seen, if any, are signals from dark matter and we use the Poisson 
statistics to find the parameter space excluded at a given confidence level 3 . Integrating 

2 There is no consensus regarding how to bin the DAMA data. The analysis of (18) for instance is based 
on the published number of events given by the collaboration (1), which are combined in two bins (2-6 keVee 
and 6-14 keVee). The analysis of (21) is based on the 36 bins given in the figure 9 of Ref. (1), recombined 
in the first 16 bins plus 1 bin from 10-20 keVee. The analysis of (20) uses all the 36 bins, while (29; 30) keep 
only the first twelve bins. Here we use all the available data, like in Ref. (20), with the reserve expressed 
in text above. 

3 More sophisticated statistical methods exist to set exclusion limits, which are especially powerful in 
the presence of few events and of an unknown background, like the maximum gap, described in Ref. (63). 
However they are not systematically adopted in the recent theoretical analyses to which we refer in this 




(3.16) 



X 2 < 65). 
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over the exposure time, the upper bounds are obtained by requiring that the total num- 
ber of events N tot is compatible with the number of observed events at 99% of CL, as in 
Ref. (30), which means that we are optimistic about the occurrence of signals for dark mat- 
ter. We have checked that our simple procedure reproduces fairly well the published 90% 
CL upper bounds in the plane Mdm — o'n plane, e.g. Refs. (20; 21; 28). In the following 
we describe the main features of the considered experiments, and specify when we use a 
different approach to obtain the exclusion limits. 

CDMS : The Cryogenic CDMS experiment at Soudan Underground Laboratory op- 
erates Ge and Si made solid-state detectors. For a heavy WIMP, the Ge data are more 
constraining: we have considered the ensemble of the three released runs, with respectively 
an exposure of 19.4 kg-day (64), 34 kg-day (65) after cuts and a total exposure of 397.8 kg- 
days before cuts from (66). The searches on Si, released in (65), are more sensitive to light 
DM candidates, with a total exposure of 65.8 kg-day before cuts. The sensitivity to nuclear 
recoils is in the energy range between 10-100 keV while the efficiencies and the energy reso- 
lution of the detectors are given in (21). Therefore in the case of the Germanium detectors 
the efficiency is parametrized as e(E R ) = 0.25 + 0.05(.Er — 10 keV)/5 keV for 10 keV < 
Er < 15 keV and e(E R ) = 0.30 for the remaining sensitivity range. The total efficiency for 
the Silicium detectors is given by e{E R ) = 0.80 x 0.95(0.10 + 0.30(£r - 5 keV)/15 keV) 
for 5 keV < E R < 20 keV and e{E R ) = 0.80 x 0.95(0.40 + 0.10(E R - 20 keV)/80 keV) 
for 20 keV < Er < 100 keV. Our conservative exclusion limits are obtained by requiring 
a total number of events N to t less than 8.4, compatible with zero observed events at 99% 
confidence level. 

XENON10 : XenonlO is a dual-phase Xenon chamber operating at LNGS. The col- 
laboration has published the data analyses in Refs. (67; 68). The total exposure is 316.4 
kg-days and 10 candidate events have been seen in the recoil energy range between 4.5-26.9 
keV. In this case we proceed with a different analysis, as in Ref. (20), using the 7 bins as 
well the bin dependent efficiencies provided by the collaboration in Table I of (67), and 
a x 2 f° r Poisson distributed data. The constraints are given by imposing the GOF be 
compatible with zero observed events at 99% CL. This experiment is sensitive to inelastic 
dark matter in a similar way as the DAM A experiment, due to the close proximity of the 
target nucleus masses. 

ZEPLIN-III : The ZEPLIN-III experiment (69) in the Palmer Underground Laboratory 
at Boulby is a two-phase Xenon chamber, so is more sensitive to heavy DM particles. The 
2008 run has a total exposure of 126.7 kg-day and has seen 7 events in the measured energy 
range (E q ) between 2-16 keVee. The quenching factor for the Xenon is parameterized as 
in (29). It is energy dependent for recoil energies up to 10 keVee, qx e = (0.142.6^ + 

paper. In the literature, Poisson statistics has been used in the analysis of Refs. (20; 30), while Ref. (29) 
and Ref. (18) use the maximum gap method. The most sophisticated analysis is that of Ref. (21), where 
they use a mixture of techniques. The results of these analyses are however consistent with each others. For 
the sake of simplicity, and because we are more focused on the features of our model than on the statistical 
analysis, we use standard Poisson statistics. We have checked that our results are consistent with those 
of the other analyses. We have also verified that using the maximum gap method affects only slightly our 
exclusion limits and that our conclusions regarding the viability of the model are essentially unchanged. 
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0.005) exp [— 0.305£'g' 564 ] and then becomes constant with a value of qx e = 0.48. The 
recoil energy Er is obtained by rescaling Er = E q /qx e - The exclusion at 99% CL is 
obtained by imposing Ntot < 16. 4 

CRESST : For the CRESST experiments at LNGS, we have derived the exclusion 
limits in both the light and the heavy WIMPs scenarios. For a light dark matter we 
have used the first data release (hereafter CRESST-I) for the prototype AI2O3 detector 
module runs (72), with a total exposure of 1.51 kg-days covering the energy range be- 
tween 0.6 - 20 keV, with 11 observed events and an energy resolution given by cr(E) = 
^(0.519 keV) 2 + (0.0408 E) 2 . The eleven events after cuts produce an upper bound on 
the total number of events of 23 at 99% CL. The collaboration has changed significantly 
the detectors for the second commissioned run of 2007 (73; 74) (CRESST-II) using CaW04 
as single crystals. The presence of the heavy nucleus of Tungsten enhances the sensitiv- 
ity to spin-independent inelastic particle scattering. The energy range of the CRESST-II 
Zora and Verena detectors is 12-100 keV, with a total exposure of 47.9 kg - day before 
cuts and an acceptance of 0.9 on Tungsten recoil. The collaboration measured 7 events, 
that provides a limit on the total rate in the modules at 99% CL compatible with N to t < 16. 

Two other experiments are potentially relevant, CoGeNT and TEXONO. The CoGeNT 
experiment (75) is based on Ge sub-keV threshold crystal detectors. The binned data 
arise from the first prototype of 0.475 kg Ge with a 8.4 kg-days exposure, in the measured 
energy range between 0.388-0.983 keVee, with a quenching factor qc e = 0.2. The TEXONO 
experiment (76; 77) runs at the Kuo-Sheng(KS) Laboratory and is again a Ge crystal with 
low threshold, 0.2 keVee. The experiment is sensitive in the energy range between 0.2 -0.8 
keVee, with the same quenching factor as CoGeNT; the detector has a total exposure of 
0.338 kg-days with an efficiency of 50%. The limits set by these two experiments are only 
relevant for a two bins analysis of the DAM A data (18-21), but to avoid the cluttering of 
the figures, we do not show them in this paper. 

4. Elastic scattering of a light WIMP 

The light WIMP scenario refers to a dark matter particle with a mass in the multi-GeV 
range. This is lighter than a typical WIMP candidate, like a standard neutralino, which 
would have a mass in the multi-10 GeV range. The possibility of explaining the DAMA 
results with a light WIMP has been first considered in (12; 13) and, for neutralino can- 
didates, in (14-16). These early proposals, which made use of the lower threshold for 
scattering of light dark matter particles on Sodium rather than, say, on Germanium, have 

4 The former ZEPLIN-II experiment is similar to ZEPLIN-III. ZEPLIN-II had a total exposure of 225 
kg-day and 29 events were observed. In the analysis of (28-30) these events are accounted as being possible 
signals of dark matter. The recent work of Ref. (70) however claims that all the events are consistent with 
being background, which obviously makes the constraint from ZEPLIN-II more stringent. In analogy with 
the Feldamn and Cousins statistics (71) they used, here we have considered a Poisson statistics compatible 
with at most 10 events which gives exclusion limits similar from those derived in (70). However we do not 
reach the same conclusions regarding the constraints on inelastic dark matter, see Section 5. 
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Figure 1: Elastic scenario The averaged modulated amplitude S m for the best fit point (see 
Figure 3). The relevant parameters are given in the table on the right. The Xmin ^ s evaluated for 
34 d.o.f. 

been excluded by more recent experimental results. However, the possible effect of chan- 
nelling in the Nal crystals of DAMA (61; 78), which effectively lowers the threshold for 
detection of nuclear recoils and reinforces the impact of collisions on Iodine, has reopened 
the possibility of explaining the data with a light WIMP (17). The impact of channelling 
has been further studied in many recent works (18-22; 62). A general conclusion of these 
works is that a light WIMP candidate may give a honest fit to the DAMA data, but is 
only marginally consistent with current exclusion limits, with only a small region in the 
mass vs (SI) cross section parameter space remaining at the 3a level. A further issue of a 
light WIMP scenario is that the total rate of nuclear recoils tends to be large and raising 
at low recoil energies, features that are absent in the data and which furthermore imply 
that background in DAMA is unexpectedly small (20; 21). Despite these caveats, there is 
much appeal to a light WIMP candidate. For instance, it is intriguing that a dark matter 
candidate with a mass of a few GeV would have an abundance similar to that of ordinary 
matter (79-85). Also, there is potentially a whole zoo of light WIMP candidates (25) which 
are yet not constrained but which are within reach of existing or forthcoming experiments. 

The elastic scenario has already been addressed in Ref. (24), following an analysis 
based on a two bins version of the recent DAMA data of Ref. (18). In Ref. (24) it was 
shown that the elastic scattering of a light scalar particle interacting dominantly through 
the Higgs channel, may be simultaneously compatible with both the DAMA data and with 
the WMAP cosmic abundance 5 . This limit of the IDM is also the simplest instance of a 
Higgs portal (87) and so has a larger scope. A singlet scalar as a candidate for dark matter 
has been discussed in various works (50-52). 

The model independent analysis of Ref. (18), which was based on DAMA data grouped 



5 This is also specific to a scalar dark matter candidate (24). For instance, it is well-known that anni- 
hilation of a fermionic singlet through the Higgs is a P-wave suppressed process. The coupling required 
to explain the relic abundance is then way too large to be compatible with direct searches (24). More 
complicated fermionic models are however viable (86). 
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Figure 2: Elastic scenario (model independent, f p — /„): Allowed regions in the plane er° — Mom 
consistent with the DAMA annual modulation signal at 90%, 99% and 99.9% CL. The solid (brown) 
curve shows the total unmodulated rate for DAMA. The escape velocity is v esc = 650 km/s. The 
dashed (black) curve is the exclusion contour at 99% CL for CDMS-Ge, while the dotted (green) 
curve is the exclusion limit at 99% for XENON10. 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

Figure 3: Elastic scenario (predictions of IDM) Left panel for v esc = 450 km/s and right panel 
for v esc — 650 km/s. Allowed regions in the plane \\h () \ — Mjj consistent with the DAMA annual 
modulation signal at 90%, 99% and 99.9% CL. The red (solid) lines delimit the WMAP 5 years 
bounds at 3 a CL (88). The solid (brown) curve shows the total unmodulated rate for DAMA. The 
dashed (black) curve is the exclusion contour at 99% CL for CDMS-Ge, while the dotted (green) 
curve is the exclusion limit at 99% for XENON10. 
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Figure 4: Elastic scenario (predictions of IDM) Allowed regions in the plane pdm— Mh consistent 
with the DAMA annual modulation signal at 90%, 99% and 99.9% CL for a fixed coupling \\h \ = 
0.12. The escape velocity is fixed to v esc — 650 km/s. The bounds are the same as in Figure 3. 

in two bins, 2 — 6 keVee and 6 — 14 keVee gave a SI cross section in the range 

3 x 1(T 41 cm 2 < erf 7 < 5 x 1(T 39 cm 2 (4.1) 

and a dark matter mass in the range 

3 GeV < M DM < 8 GeV, (4.2) 

where the upper bound on the mass is set by the exclusion limit from XENON. Subse- 
quently, various groups have emphasized the relevance of using the full set of spectral data 
provided by DAMA (19-22). These spectral data are shown in Figure 1, together with the 
theoretical spectrum for the elastic scattering of a 12 GeV candidate. In the IDM, this 
would correspond to a light scalar with coupling \\h \ = 0.24 to the Higgs, and = 120 
GeV. The raise at small recoil energies is typical of elastic scattering of a light WIMP. 

When the full spectral information of DAMA is used, the allowed region of parame- 
ters is much reduced compared to the ranges of Eqs. (4.1) and (4.2). Our main results 
regarding a light WIMP in the IDM are summarized in Figures 2 and 3 . Figure 2 is model 
independent, the only hypothesis being that the coupling to nuclei may be approximated 
with f n = f p . It shows the region allowed by DAMA, up to 99.9% CL, together with the 
exclusion limits from XENON-10 and CDMS-Ge, which are the only two experiments rele- 
vant in these ranges. Our results are consistent with those of Refs. (19-22). A small region 
of masses and cross sections is compatible with all experiments, albeit only at the 3a level. 
Figure 3 gives the allowed parameter space for the IDM, or equivalently for a light scalar 
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singlet, the left panel is for an escape velocity of 450 km/s and the right one for v esc = 650 
km/s. The mass of the Higgs has been fixed to Mh = 120 GeV, but all the predictions 
only depend on the ratio A^- o /M^ and may thus be trivially rescaled (24). Also shown in 
the figure is the predictions for the abundance from thermal freeze-out, compared with the 
WMAP data, 0.091 < ^da//i 2 < 0.129 at 3a (19). We should perhaps emphasize the fact 
that no further adjustment of parameters is done in order to fix the relic abundance. 

The overlap between the region allowed by DAMA and the WMAP abundance is 
amusing, but should perhaps not be taken too seriously, independently of the fact that this 
region is actually excluded by the other experiments. One may envision at least three ways 
in which the regions of parameters may be modified. First the local abundance is fixed 
here to the fiducial value pdm = 0.3 GeV/cm 3 . Decreasing the local abundance will move 
the DAMA region and the excluded region to larger quartic couplings, without affecting 
the cosmic abundance. See also Figure 4, where we show this variation for a fixed quartic 
coupling. Second, one should take into account the fact that the coupling of the Higgs to 
a nucleon N, or ghNN, is not well-known. This is parameterized by a form factor / using 
the trace anomaly, fmj^ = (N\ J2q m q'j'l\N) = ghNNV- The reference value for / = 0.3 
(89-92), but it may vary within a quite large range 0.15 < / < 0.6 (24). Changing the 
parameter / will move around the DAMA and excluded region, without affecting the relic 
abundance. Finally, one may modify the relic abundance by changing the cosmology. 

5. Heavier WIMPs and inelastic scattering 

In the inelastic Dark Matter scenario, or iDM for short, an incoming dark matter particle 
DM\ from the halo is supposed to scatter with a nucleus into a slightly heavier state DM2, 
with a mass splitting 5 = Mdm 2 — M-dm x ~ 100 keV. In this scenario, which has been first 
proposed in Ref. (26) and confronted to the recent data in Refs. (27-31), a much broader 
range of dark matter candidates may both fit DAMA and be consistent with the other 
experiments. To the naked eye the fit to the modulated spectrum is better than in the 
elastic case, see Figure 5, although the GOF to data is comparable. On one hand, the iDM 
scenario opens the exciting possibility that there might be in the dark sector more than 
just one plain particle. On the other hand, from the point of view of particle physics, the 
naturalness of the small mass splitting needs to be addressed (30). 

To illustrate the relevance of the inelastic scenario, we show in Figure 6 four snapshots 
of the Mdm — o"n plane for different mass splittings, using Eq.(3.4) with f p = f n with 
increasing splittings. These panels show that, as the mass splitting is increased the two 
regions of Figure 3 (which correspond to zero mass splitting) merge and move toward higher 
dark matter masses and larger cross sections, and more solutions open. 

The possibility of dark matter from a weak SU{2) doublet has already been considered 
in the literature (27; 28; 30). Here we focus on the prediction of the IDM, which is the only 
model among Scalar Multiplet Models where a mass splitting of the order of 100 keV, small 
enough to be relevant for inelastic scattering, can appear at the renormalizable level. Notice 
that when 5 = Ma q — Mh *C Mh , the only accelerator constraint on IDM comes from the 
Z boson decay width. In particular, this limit is not excluded by LEPII measurements (48). 
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Figure 5: Inelastic scenario (predictions of IDM) On the left, the averaged modulated amplitude 
S m for the best fit points (for v esc = 600 km/s, see also Figure 8). On the right: Values for the 
global and local minima of the chi-square distribution for the best fit points. The Xmin ^ s evaluated 
for 34 d.o.f. . 

In the IDM, an accidental U{\)pq symmetry arises in the limit of vanishing A5, so that 

the mass splitting 8 is protected against radiative corrections. In this sense, the option of 

having a small quartic coupling A5 is technically natural. 6 

For the sake of comparison with other works, in particular Refs. (29; 30) with which we 

agree, we show in Figure 7 the allowed regions in the 5 — o~ plane, where is the WIMP 

cross section with a neutron for SU (2)-type couplings (f p = 4 sin 2 6w — 1, and /„ = 1) but 

for an arbitrary scale (the horizontal line at ~ 7 • 1CF 39 is the SM Z-boson exchange) 
7 

In Figure 8, we show the predictions of the IDM, the left panel is for an escape velocity 
of 450 km/s and the right one for v esc = 600 km/s. These figures show that there exists 
a whole range of candidates, between Mh ~ 535 GeV and M# ~ 20 TeV, which are 
compatible with DAM A and all the other experiments, and which have a relic abundance 
consistent with WMAP. The candidates in the pink region (light grey bottom) have an 

6 As put forward in Refs. (44; 45), the Inert Doublet Model may be embedded in SO(10), provided 
the Standard Model fermions are also odd under Z2 or matter-parity, with the Standard Model in a 10 of 
SO(10) and the inert doublet in a 16. A very interesting feature of this framework is that the coupling A5 
is absent at the level of renormalizable operators (45). Higher order operators may break the accidental 
PQ symmetry while preserving matter-parity (which is a gauged discrete symmetry and thus protected). 
In Ref. (45), breaking by quantum gravity effect are estimated to give A5 ~ (M so (10) /Mpi) n with n > 1. 
For n = 2 and Mso(io) ~ 10 16 GeV for the scale of grand unification scale, gives A 5 ~ 10~ 7 , in the range 
required to explain DAMA through inelastic scattering, Eq.(2.4). We would like to emphasise that higher 
operators may also arise within SO(10). For instance, with scalar fields in a 126, which is required to 
break B — L, one may build a dimension 5 operator suppressed by some scale A (at tree level, this may be 
mediated by scalars in a 54), which gives A5 ~ U126/A after breaking of B — L (we thank Thomas Hambye 
for discussions on this point). 

7 Unlike Ref. (70), we find that ZEPLIN-II is not much more constraining than other experiments, like 
ZEPLIN-III, in the region of the inelastic scenario. This is due to the fact that the DAMA region we have 
found, which is consistent with results in Refs. (29; 30), differs from that of Ref. (70). 
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Figure 6: Inelastic scenario (model independent, f p — f n ) Allowed regions in the plane Mbm— &n 
consistent with the DAMA annual modulation signal at 90%, 99% and 99.9% CL for different mass 
splittings. From left to right and from top to bottom, respectively, <5 = 5 keV, 8 = 20 keV, 8 = 80 
keV and 8 = 150 keV. The escape velocity is v esc — 650 km/s. The total unmodulated rate expected 
at DAMA is given by the solid (brown) line. The dotted (green) curve is the exclusion contour at 
99% CL for XenonlO, the dashed (black) curve is the upper bound coming from CDMS-Ge, while 
the dotted-dashed (blue) curve is the exclusion limit at 99% for CRESST-II. 

abundance below the WMAP observation. The grey region (top) corresponds to the Griest- 
Kamionkowski unitarity limit Mh < 58 TeV, see Refs. (43; 93). Notice that the cross 
section in the inelastic dark matter scenario is fixed, as the scattering occurs dominantly 
through a Z boson. Therefore the relevant parameter space is Mh — 5. The brown (dark 
grey bottom) region is excluded by LEPI measurement of the width of the Z. 

Heavy candidates with a mass above 535 GeV have been shown to be compatible 
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Figure 7: Inelastic scenario (model independent, f p — 4 sin 2 9w — 1, f n = 1) Allowed regions in 
the plane S —- consistent with the DAMA annual modulation signal at 90%, 99% and 99.9% 
CL for a fixed mass Mum — 525 GeV. The escape velocity is v esc = 600 km/s. The dashed 
(black) curve is the exclusion contour at 99% CL for CDMS-Ge, the dotted-dashed (blue) curve is 
the exclusion limit at 99% for CRESST-II, the long-dashed (red) line denotes the exclusion limit 
coming from ZEPLIN-III and the short-dashed (red) line is that from ZEPLIN-II. The grey solid 
line corresponds to the Z exchange with a neutron. 

with WMAP in Ref. (43), irrespective of the precise value of the small splitting 5. The 
agreement can be achieved by tuning the values of \X2 and M H + , which control the amount 
of coannihilation. For candidates with a mass lower than 535 GeV, direct annihilation 
channels into Z or W bosons, when opened, reduce the relic abundance below the WMAP 
value. For candidates lighter than the W threshold, but above the lower bound set by the 
LEP measurement on the decay width of the Z, the coannihilation cross section through 
the Z boson in the s-channel is large due to the proximity of the Z pole. Even at Mjj = 
80 GeV, the total coannihilation cross section is about 25 pb, much larger than 1 pb, the 
typical value needed for WMAP. Therefore, middle mass range candidates with inelastic 
interactions that are consistent with DAMA have a small relic abundance. 

With respect to this problem, we put forward two options. One possibility is that the 
neutral scalar is a subdominant component of dark matter and we show how the fit changes 
as the local density of the neutral scalar is changed 8 in Figure 9. Alternatively one may 
envision the possibility that there is a charge asymmetry in the dark sector. Concretely, 
we may use the approximate PQ symmetry, which becomes exact in the limit A5 — > 0. 

8 The local density can also vary due to the clumpy nature of DM in galactic halos. The value of 
0.3 GeV/cm 3 usually quoted only reflects an average density of the smooth component at the location of 
the Earth. 
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Figure 8: Inelastic scenario (predictions of IDM) Left panel for v esc — 450 km/s and right panel 
for v esc = 600 km/s. Allowed regions in the plane <5 — Mh consistent with the DAMA annual 
modulation signal at 90%, 99% and 99.9% CL. The solid (brown) curve shows the total unmodulated 
rate for DAMA. The dashed (black) curve is the exclusion contour at 99% CL for CDMS-Ge, while 
the dotted-dashed (blue) curve is the exclusion limit at 99% for CRESST-II. The brown (dark grey 
bottom) region is excluded by LEPI measurement of the width of the Z; the pink region (light 
grey bottom) is below the WMAP abundance; the grey region (top) corresponds to the Griest- 
Kamionkowski unitarity limit, Refs (43; 93). 



We would like to argue here that, for small enough A5, an initial PQ charge asymmetry 
produced in the early universe may survive until after freeze-out. Indeed, in the limit 
A5 — > 0, the complex, neutral scalar field H n = (Hq + iAo)/V%, has a PQ charge +1 
which is conserved by gauge interactions and the scalar interactions given by the potential 
Eq. (2.1). For small A5, this charge is only broken by processes proportional to A§, in the 
annihilation channels H n H n — > hh or H n H n — > {h) — > //, and the corresponding scattering 
processes H n X — > H*X. In full generality, the resolution of the Boltzmann equation that 
governs the evolution of the comoving charge density A = (n# n — nu*)/s(T), where s(T) 
is the entropy density at temperature T, are likely to be intricate. Details about the 
electroweak phase transition, like the temperature dependence of the Higgs mass or the 
fate of the Goldstone bosons, can affect the result. Also, the number of processes that 
violate the PQ charge is temperature dependent. In particular, the annihilation channels 
into Higgs, top quark and bottom quark may or may not be present. 

To simplify the problem, here we will consider that the contact interactions H n H n — ^> 
hh (annihilation) and H n h H*h (scattering) are the dominant processes. This assump- 
tion is justified by the fact that channels with an intermediate Higgs like H n H n — > {h) — > ff 
are further suppressed by Yukawa couplings. Also, at high temperature before the freeze- 
out, njj n — nu*. With these hypotheses, the evolution of the asymmetry is controlled by 
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Figure 9: Inelastic scenario (predictions of IDM) Allowed regions in the plane Pdm ~ Mh con- 
sistent with the DAMA annual modulation signal at 90%, 99% and 99.9% CL for a fixed mass 
splitting S — 120 keV. The solid (brown) curve shows the total unmodulated rate for DAMA. The 
escape velocity is v esc — 600 km/s. The dashed (black) curve is the exclusion contour at 99% CL 
for CDMS-Ge, while the dotted-dashed (blue) curve is the exclusion limit at 99% for CRESST-II. 



a simple Boltzmann equation 

dA ( Y n Y h \ . . . 

with a = x/Hs(T), x = Mdm/T, H is the Hubble parameter, and Y{ = Hi/s(T) is the 
comoving number density of the species i. The thermal average of a cross-section a for the 
process 12 «-> 34 is given by (av) = 7(12 ^ 34)/nf (T)n 2 q (T), with (94) 

7 (12 <-» 34) = J ds yfs Kti^s/T) p 2 lcm a(12 <-» 34) . (5.2) 

For Mom 0, the equilibrium number density reduces to 

n^ = 9A $f r K 2 (M DM /T)^ 9 - 1 ^ , (5.3) 

with g the number of internal degrees of freedom. 

To further simplify, we will neglect the masses of the particles in computing thermal 
averaged quantities. This will provide an upper bound on the relevant interaction rate, as 
a finite mass can only reduce the number density in the thermal bath. In this limit, the 
annihilation and the scattering cross sections are simply given by 

aA = 1 r s = S- s ' (5 - 4) 
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with s the center of mass energy and the thermal average 7 is given by 



1A = 




(5.5) 



2 8 7T 5 ' 



Using the relation H = 1. 66 ^/g^T/Mpi, the out-of-equilibrium condition for the interaction 
rate, Ta = {(JAv)n eq < H, implies 



By comparing this relation with Eq. (2.4), we see that PQ breaking processes are indeed 
out-of-equilibrium in the early universe, for candidates lighter than about 100 GeV. For 
heavier candidates [Mom — 1 • • • 10 TeV), the charge asymmetry may be washed out, a 
thermal relic abundance consistent with WMAP can be obtained through the standard 
freeze-out mechanism, without invoking an asymmetry. 

6. Conclusions and Prospects 

In this article we have confronted the inert doublet model or IDM to DAMA. In this very 
simple framework, which consists of just one extra, inert, Higgs doublet, it is possible to 
explain the DAMA data in two particular limits of the model. 

In the first limit, which is protected by a custodial symmetry, all, but the lightest, 
scalars are decoupled and the model is essentially equivalent to a singlet scalar extension of 
the Standard Model in which the dark sector interact with the rest of the world through the 
Higgs channel, aka the Higgs portal. A fit to the DAMA data requires a rather light scalar 
dark matter candidate, in the multi-GeV range, which undergoes elastic, spin independent 
scattering with nuclei (essentially Iodine) in the detector. This scenario, which works if 
channelling is effective, is severely constrained by exclusion limits imposed by XENON-10 
and CDMS. There is also some tension with the fact that the total, unmodulated rate needs 
to be quite large, possibly saturating the signal in the detector. The analysis presented 
here is essentially an update of Ref. (24), were it was shown that, without adjusting any 
parameter, such a model may simultaneously explain the DAMA data and be consistent 
with WMAP. Our more refined analysis, which is based on using the full spectrum data 
provided by DAMA, concurs, as is shown in Figure 3. Although this scenario, which 
encompasses a potentially broad range of models, is challenged by the other experiments, 
we would like to re-emphasize here that it predicts numerous other signatures. This stems 
both from the larger abundance of a light WIMP, and from the rather large couplings or 
cross section required to fit the DAMA data. Depending on the abundance at the Galactic 
center, one may have a rather large flux of gamma rays from dark matter annihilations, 
which is already constrained by EGRET, and should be within the reach of the forthcoming 
Fermi/GLAST data (24; 25). Capture by the Sun and the subsequent annihilation into 
neutrinos may be constrained by Super-Kamiokande data (21; 22; 39; 95). Annihilations 
in the Galaxy may produce antimatter, in particular in the form of anti-protons and anti- 
deuterons (40; 62). Finally, a light WIMP in the form of a scalar coupled to the Higgs 




(5.6) 
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would imply the Higgs mostly decays into a pair of dark matter particles, with striking 
consequences for its search at the LHC (24; 51; 52). It is also intriguing (albeit puzzling) 
that in this scenario the abundance of dark matter is similar to that of baryons. 

In the second limit of the IDM considered in the present paper, the two neutral scalars 
of the inert doublet are almost degenerate, with a splitting 0(100 keV), necessary to 
explain the DAMA data through the inelastic scattering of the lightest of the partners, the 
so-called inelastic Dark Matter or iDM scenario. This limit of the IDM is also protected 
by a symmetry, a U{\)pq in this case. This implies that the splitting, although fine tuned 
and requiring a small parameter in the form of a quartic coupling, is technically natural. 
Here we show that this very simple model may explain the DAMA data and be consistent 
with WMAP data for a whole range of candidates, with a mass between ~ 535 GeV and 
~ 50 TeV. This is shown in Figure 8, in the mass/mass splitting plane. This scenario 
has a GOF which is comparable to that of the elastic scenario, but it is less constrained 
by the other experiments even if, as far as we know, it does not have a wealth of other 
potential signatures, expected possibly neutrinos from the Sun (96; 97), or using specific 
direct detection setups (31). An issue with the iDM scenario is that it is not easy to 
explain the smallness of the mass splitting required to fit the data (26; 30). In the IDM, 
there is no fundamental explanation for this smallness, but higher scalar multiplets may 
do better (43). For instance, an inert scalar SU{2) triplet with hypercharge Y = 2 and 
with coupling only to the SM Higgs, has two neutral scalars which are degenerate at the 
level of renormalizable operators. A splitting may arise through a dimension 6 operator, 
and thus is suppressed by a factor ~ v 2 /A 2 , with A>ua new scale. Another interesting 
possibility opened by the model discussed here is that the small mass splitting may be 
related to a tiny violation of a global U{\)pq symmetry. This implies that, at least for 
not too heavy candidates, Mqm ~ 100 GeV, an initial U{\)pq asymmetry in the Early 
Universe may survive and not be washed-out before annihilations of dark matter particles 
actually freeze-out. An asymmetry in the dark scalar sector could for instance be generated 
through leptogenesis, if right-handed neutrinos, odd under Z2 are also introduced in the 
model. The addition of right-handed neutrinos is actually a very natural extension of the 
IDM (34), which might also explain the origin of the SM neutrino masses as being due 
to radiative corrections. There is some tension between the fact that leptogenesis requires 
rather heavy right-handed neutrinos, Mjy ~ 10 11 GeV, while in the model we have in mind, 
the SM neutrino masses are expected to be 

5 • M DM 

with 5 ~ 100 keV and Mdm ~ 100 GeV, which a priori gives a too small prediction for 
m v . This holds if we assume that the heavy right-handed neutrinos have similar masses 
and similar Yukawa couplings to SM neutrinos. As emphasized in (43) this does not have 
to be the case, and a hierarchical structure of Yukawa couplings together with a light 
right-handed neutrino, Mjv ~ 1 TeV, might lead both to successful leptogenesis and SM 
neutrino masses in agreement with observations. We leave this interesting possibility for 
future investigations.. 
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